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We examine the subtle differences between two possible frameworks for the calculation of quarko- 
nium production. The differences between the two methods are not a concern when one calculates 
only to leading order in the relativistic expansion, but when relativistic corrections are included, the 
two formulations seem at first to be at variance with one another. They can be reconciled however 
via a relativistically corrected mass relation relating the boundstate mass and the quark mass. 



I. INTRODUCTION 



We present here an observation regarding relativistic corrections in factorization formulas for quarkonium produc- 
tion, as calculated in the non-relativistic QCD framework JlJ . We point out that one can conceive of two subtlely 
different frameworks in which calculations of quarkonium production might be carried out. We illustrate the algorithm 
for both approaches through the detailed presentation of a simple example. We explain that the results of the methods 
differ on two points: on the question of the definitions of the NRQCD matrix elements; and on the question of the 
mass parameters that appear in NRQCD factorization formulas. In Method I, factorization formulas are computed in 
such a manner that only the quark mass m c appears, while in Method II, the formulation generally entails the mass 
of the quarkonium boundstate (which we will denote with the upper case Mh) as well as m c . 

The difference between the two methods is not a concern when one calculates quarkonium production and decay 
rates only to leading order in the relativistic expansion. In that approximation, one sets Mh = 2m c , and the two 
methods give compatible results. However, when relativistic corrections are included (and when Mh no longer is 
equatable to 2m c ), we find that the two formulations seem at first to be at variance with one another. 

The results of the two methods can be reconciled using a relativistically corrected mass relation, which is an equation 
relating m c , Mh and the NRQCD production matrix elements. 

A variant of the relativistically corrected mass relation has recently been discussed in Ref. [||. The appearance of 
this work increases the possible interest of the NRQCD community in our observation. 

In order to illustrate our observation, we take the simplest possible example, that of charmonium production from 
the collision of two a ficticious elementary (colorless) scalars, denoted by <f>: 

<f> + cf> — > charmonium. 



The underlying short-distance physics is the partonic process 

+ <t> — ► C + C 

which is described by the interaction term 

dnt = g4>{x)<l){x)li{x)^¥{x) (l.l) 

where ty\x) is the four-spinor field operator for the charm quark (with color index i), and where g is a coupling 
constant of mass dimension —1. The reaction <f> + (f> — > c + c produces a cc pair in a color-singlet state with angular 
momentum 2S+1 Lj = 1 Sq. This cc can hadronize into the r\ c charmonium particle. We choose this process because 
of the simplicity of the spin structure and because there are no color-octet contributions which would distract from 
our main point. 

The determination of production factorization formulas in NRQCD is always based, in one way or another, on the 
Feynman amplitude M for the production of a free heavy quark and free heavy antiquark. In particular, we desire to 
know how this quantity depends on q, the relative three momentum of the quark and antiquark in the quark-antiquark 
restframe. 

In our toy example, the Lorentz-invariant Feynman amplitude is given by 

M(cc(d)) = guv = 2gE q £Cr- (1.2) 
(Colour indices have been suppressed.) The Dirac four-spinors u and v are normalized so that uu — vv — 2m c . 



The passage from the middle piece of Eq. 1.2 to the right-hand side is an exact relation, and not a nonrelativistic 



approximation. It entails recasting the four-spinors (from the Pauli-Dirac representation) in terms of the two-spinors 
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£ and £, as described in Ref. M. £ is the two-spinor for the heavy quark, and £ (called r) in Ref. S) is the two-spinor 
for the heavy antiquark. a and r label heavy quark spins; they are not spinor component indices. The two-spinors 
are normalized so that £j£ r = 6 aT . As to E q , this is simply a short form for \/ q 2 + m 2 . The invariant mass of the 
system consisting of the pair of free on-s hcll heavy quarks is 2E g . 



The spinor-structure given in Eq. 1.2 is such that multiplying the expression by the Clebsch-Gordon coefficients 



(1/2, a; 1/2, r|0, 0) and then summing over spins gives a non-zero result. Howe ver, doing the same with the Clebsch 



Gordon coefficients (1/2, a; l/2,r|l,M) gives zero. We therefore see that Eq. 1.2 represents the production of a cc 
pair in a spin state J,Mj = 0, 0, as was previously stated. 

II. RELATIVISTICALLY CORRECTED FACTORIZATION FORMULAS: METHOD I 

A. The Concept of Matching 

The NRQCD factorization formalism exploits the equivalence of perturbative QCD and perturbative NRQCD. The 
computation of the factorization formula for the production of a quarkonium boundstate starts with a computation of 
the production rate for free on-shcll heavy quarks and antiquarks; this computation is performed in some perturbative 
model of underlying short-distance physics, which could be perturbative QCD, electroweak physics, or (as here) a toy 
model. In a separate calculation, one uses NRQCD perturbation theory to compute the free-quark matrix elements 
of the four-fermion operators of the NRQCD effective theory. The perturbative equivalence of NRQCD and the full 
underlying theory is expressed by the matching condition (Eq. (6.7) of Ref. [[l]) 

-(-(q))L ort -E^^(°l^ ) l°)|pert NRQCD . (2-1) 

pel l TTlc 

The left-hand side is to be computed in the perturbative model of underlying physics. It is the rate for the production 
of a heavy quark pair, written expressly as a function of q, the relative three-momentum of the c and c in the cc 
rcstframe. As well as being a function of q, the left-hand side is also a function of whatever other kinematic variables 
are required for parametrizing the process at hand. As to the right-hand side, it is simply a linear combination 
NRQCD free-quark matrix elements (O|0n C< " <? ' > |0) such as those defined in Appendix I. The sum over n in the right- 
hand side is a Taylor series; each of the matrix elements is a function of q and m c only, with the power in q 2 /to 2 , 
increasing with each higher term. The matching condition is intended to allow a determination of the short-distance 
coefficients F n . By definition, these are independent of q. How the matching condition is to be employed will become 
clear in the example below. The basic idea is that, by calculating the production rate in the full underlying theory 
on the left-hand side, and by massaging it into a Taylor series of the form of the right-hand side, one can then read 
off the F n . d n is the mass dimension of the product of the field operators in O n . 

Once the short-distance coefficients F n have been calculated for free on-shell heavy quarks in the perturbative 
theory, the rate for the production of some quarkonium boundstate H is then known to be (Eq. (6.4) of Ref. jlj]) 

^) = E^Wl°)- ( 2 - 2 ) 

The above expression is a factorization formula. The NRQCD non-perturbative matrix elements (0|O^|0) are to be 
thought of as empirical parameters. A fundamental feature of the NRQCD factorization formalism is that the coeffi- 
cients F n are universal to both free quark and boundstate production. In general, the short-distance coefficients and 
the hadronic matrix elements also depend on the NRQCD cut-off, but we will not bother to display this dependence, 
since it is not a point of focus in our discussion. 

B. Manipulation of Phase Space for Method I 

In Method I, we interpret the matching condition, Eq. |]l], to involve (on both sides of the equation) the rate of 
production of free heavy quarks; this interpretation defines Method I. In such an approach, the physical meaning of 
<r(cc(q)) is fixed via the concept of the total cross-section for free on-shell cc production: 

-i'-'-S = / 4% <Kcc(q)) • (2.3) 
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By construction, <r(cc(q)) depends on q, m c and other kinematic variables of the free cc production process, but it 
logically does not depend on the mass of quarkonium boundstate, Mjj- 

Our goal in this section is to obtain a useful general expression for cncc(q)) for the purposes of matching in Method 
I. We procede by first computing the production rate for a pair of free on-shell quarks. The total rate is given by the 
well-known formula 



(2.4) 



where, since we are only considering spin-summed rates, we sum over all quantum numbers of the final state free 
heavy quark and free heavy antiquark. Here, X represents the sum of the momenta of all final state particles other 
than the two heavy quarks. ^2 X represents the measure of phase space integration over all such momenta. 
Let us now define the total and relative four- momenta of the heavy quarks: 



P = Pc 

2Q = Pc 



Pc 
Pc- 



P is the total four-momentum of the cc system in the lab frame. Q is the relative four-momentum of the c and c in 
the lab frame. The lab-frame four-momentum Q is related to the rest-frame three-momentum q by the relation 

Q M = AV 

where is the Lorentz boost matrix connecting the lab frame with the cc restframe. Recalling that the invariant 
mass of the cc system is 2E q = + q 2 , we see that the components of the total four-momentum of the cc system 

are 



P = {y/4E% + P 2 ,p} 



We now manipulate the phase space integral in Eq. 2.4 by defining the new variables of integration 



P = p c 
2Q = p c 



Pc 
Pc 



With these new variables, the total rate is written 



o ( cc ) 



The Jacobian for the above variable redefinition is unity. 

We now change variables of integration once more, replacing Q % in favor of q 3 . Concerning this transformation, we 
have the following exact identity: 



d 3 Q 
d 3 q 



1 



2p° 2 V l 
dQ l 



Dot 



dq j 



1 



2E n 



1 - 



2p° 2p a - 
P 2 q 2 

i2(p°) 2 £; 2 



(P 



0^2 



2„2 



1 - 



P 2 q 



12(P ) 2 E 2 



so that the total rate can be written in the form: 

d 3 q d 3 P 



<j[cc) 



(2tt) 3 {2ir) 3 2E q P° ^ flux 



E W- Si (P> -X-P) (2^) 4 \M{cc( q ))\ 2 



We are now able to extract a(cc(q)) , the rate as a function of q. This is seen to be 

- J (2^7^ 5 i ^ - x - p ) ^ i^wy i 2 



(2.5) 
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C. Matching Condition for Method I 



Eq. 2.5 gives an expression for the left-hand side of the basic matching condition, Eq. 2.1. The matching condition 



now becomes 

/ j^S^ £ as ^ -x-n w i^PMi 2 = ^ WW • (2.6) 

We will take this as our matching condition for Method I. It is to be used as follows: Keeping in mind that E q and 
P° are functions of q, one computes the left-hand side in the underlying perturbative theory; one expresses this as a 
Taylor series in q 2 /m^, and then massages it into the form of the right-hand side; one then reads off the F n . Because 
of the universality of the F n , the F n determine d in the matching procedure also serve in the factorization formula for 
the production of boundstate quarkonium, Eq. |2.2| , which is 

In the next section, we present an illustrative example of the method. 



D. Example of Calculation in Method I 



We now consider our toy example, </>+ (j> — * r\ c . The underlying interaction is expressed in Eq. 1.1. Let us compute 
the left-hand side of the matching condition Eq. 2.6. In this example, there is no sum ^2 X - We have 



S 3 (P)S(E f - Ek)(2*)* |-M( C c(q))| 5 



( 27 0^I 8{E f - Ei) |A<(cc(q))| s 



(2^)^j 5(E f - Ei) 4E 2 q g 2 N c £ 



^g 

Am 2 



S(E f - Ek) 



mi 



(2.7) 



where we have used the Feynman amplitude given in Eq. 1.2. This can be written in terms of the free-quark NRQCD 
matrix elements given in Appendix I: 



lhs 



4m? 



5(E f - E^ 



f^Cs )\o)-^(o\vr w Cs )\o) + 



1 



(2. 



We infer from this that the production rate for the charmonium boundstate r\ is given by 



a(r)) 



Kg 

4m; 



S(E f - Ei) 



(0\OieS )\0)-—(0\V?eS )\0) 



(2.9) 



This is the factorization formula for Method I. It corresponds to Eq. (6.8a) of Rcf. fjj. We have included the leading 
piece and the first relativistic correction. 

We must reemphasize that the factorization formulas in this approach make no explicit reference to the boundstate 
mass Mh- The short-distance coefficients can logically depend only on m c . 
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III. RELATIVISTICALLY CORRECTED FACTORIZATION FORMULAS: METHOD II 



A. An "Effective Amplitude Squared" for Quarkonium Production 

Another approach to carrying out calculations of quarkonium production rates is suggested in Ref. S. It is 
our interpretation that in this work, it is implicitly assumed that there exists a meaningful theoretical quantity - 
an "effective Feynman amplitude (squared), " \M(H)\ 2 — for the production of a quarkonium boundstate. This 
assumption is tacitly present in the first equation in Ref. , which gives the cross-section for quarkonium boundstate 
production in terms of \M(H)\ 2 : 

Here, X represents the sum of the momenta of all final state particles other than the quarkonium boundstate. 
represents the measure of phase space integration over all such momenta. It is important to appreciate that here Pjj 
represents the four-momentum of the quarkonium boundstate. The components of Ph are {Eh,Ph}-, where Eh is 
defined as yf M H + P H and where Mh is the mass of the boundstate. 

Ref. H postulates that the production rate for quarkonia can be written in a factorized form 

We use the symbols (( and )) for the bras and kets of Method II as a reminder that here the matrix elements ((0|O^|0)) 
are defined differently from the (0|O^ |0) of Method II. The reader must examine Appendices I and II for the details 
concerning the different conventions and definitions. The distinction concerns the normalization of states. The 
Fn(?n c , Ph) are the short-distance coefficients of Method II, distinct from the F n of Method I. With the states being 
normalized differently, the NRQCD matrix elements (and also the short-distance coefficients) arc of different mass 
dimension in the two methods. 

B. Matching Condition for Method II 

In order to determine the short-distance coefficients T n for a given quarkonium production process, we require a 



matching condition, analagous to Eq. 2.6. The condition proposed in Refs. |3| and m is 



£ (2tt) 4 5\ Pl -X- P(q)) \M(cc( q ))\ 2 = £ (0|O* '>|0> (3.3) 

X n m c 

This matching condition lends itself to being interpreted in the following manner: one computes |Al(cc(q))| 2 in the 
underlying perturbative theory; this allows one to express the complete left-hand side as a Taylor series in q 2 /m 2 (it 
must be kept in mind that P° depends on q); one then massages the resulting expression into the form of the right- 
hand side; finally one reads off the T n . These T n serve in the factorization formula for the production of boundstate 



quarkonium, Eq. 3.2 



In the next section, we present an illustrative example of the method. 

C. Example of Calculation in Method II 

We now calculate the factorization formula for our toy example using Method II. The left-hand side of the matching 



condition, Eq. 3.3, is 



lhs = 5 3 (P)8(E f - Ek) (2tt) 4 \M{cc(q))\ 2 

= 6 3 (P)6(E f - Ek) (2tt) 4 4g 2 E 2 q N c £ £tr (3.4) 

This can be written in terms of the free-quark NRQCD matrix elements given in Appendix II: 
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lhs = g 2 6 3 (P)6(E f - Ei) {2irf (0\of g) C S )\0) (3.5) 

Interestingly, the Taylor expansion stops after the first term, but this is not a generic fea ture of this method. In 
general, there would be an infinite series of terms in powers of q 2 /m^. We infer from Eq. |3 . 5| that the production rate 
for the charmonium boundstate rj c is given by 

°(n) = J ^^Eh 9 2 S 3 (Pu)S(E f - Ei) (2^) 4 (O|0?(^o)|O» 

:T,/2 {0\&ieS )\0}6(E f -E t ) (3.6) 



2M* 



This is the factorization formula for Method II. 

IV. COMPARISON OF THE TWO METHODS 

It is instructive to gather together the matching conditions for the two methods, for the purposes of comparison. 

E - x - p (q)) ^- ( 2 -) 4 i-^(-(q))i 2 - E »f q) \o) (4.1) 

where in both cases, the four-momentum P(q) is a function of q in the sense that, while P is strictly independent 
of q, P°(q) is given by P°(q) = \J AE q + P 2 . This last point is not an issue in the toy model, but is an issue for 
the more general situation in which (the total initial energy) is fixed (and is independent of q) , as in the case of 
b -> J/ip + s. 

Once the short-distance coefficients have been found with these matching conditions, the total production rate is 
given by 

^) = E^t<°i^i°> 

=j: rn{ 7^ H) mom (4.2) 

The reader will appreciate, by trying a few simple examples, that the state normalizations and conventions in 
Appendix I are ideally suited to Method I, and cannot be naturally employed in Method II. The state normalizations 
and conventions in Appendix II are ideally suited to Method II, and cannot be naturally employed in Method I. 

V. RECONCILIATION OF THE TWO METHODS ? 

We have obtained the following factorization formulas for the charmonium production process <j) + ^ — * V' 



2 

Method I: a(n)= ^ 5(E f - EA 
Ami 



(0|^( 1 ^ )|0)-^(0|^( 1 5o)|0) 



Method II: a( V )= ^ 6(E f (0\OJC S )\0) (5.1) 
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The reader can verify that for charmonium production from the decay of the 6-quark, b — > s + J/ip, one obtains (for 
the color-singlet contriubtions only ) 



G 



K m 3 

Method I: ct(JM) = — 
w v ' (24tt) m c 



Method II: a(J/ip) 



K ml 



(967r) 



(0 ;/ *fSi))-^(^fSi)> + ' 

v 6 TO~ 

(0^(35 1 )»- § ^(^( 3 5 1 )»+- 



(5.2) 
(5.3) 



where -ft" is some constant of mass dimension —4 involving G 2 F and CKM matrix elements. (We have simplified things 
by assuming, wrongly, that m c -C mj.) 

We have checked by explicit calculation that if relativistically corrected factorization formulas for color-singlet S- 
wave production are calculated for any quarkonium production processes, in both methods, one finds that the results 
can always be reconciled, up to relative order w 4 , using the formulas 



and 



M H = 2m c 1 



2m 2 c (0 H ) 



2M h {O h ) = (O h ) 
2M H (V H ) = (V H ) 



(l + 0(v 4 )) =2m c 1 + 



{P H )) 




(5.4) 
(5.5) 



(5.6) 



where v is the characteristic relative velocity of the heavy quarks in the boundstate. 

The first two equations seems reasonable, since one might think that they simply expresses the difference between 
the nonrelativistic state normalizations of Appenidix I and the relativistic state normalizations of Appendix II. As to 
the second equation, it also seems reasonable, since it expresses the idea that relativistic corrections to the mass of 
the boundstate are of relative order 

(V*) , 



\{o H ) 



The relativistically corrected mass relation, Eq. 5.6, trades relativistic corrections to the mass of the quarkonium 
boundstate for relativistic corrections involving the four-fermion production matrix element (V H ). Such an idea has 
recently been presented in Ref. 0. These authors use the equations of motion of the heavy quark field operators 



D 2 



to re-express the matrix element (7- >7? ) so as to obtain 

(Ol^CsyiO) = m c {M v - 2m c ) (0\CF C S a )\0) 

Solving for M v , we find 



(5.7) 



(5.8) 



Mr, = 2m r 1 



2m 2 c {O r i) 



which is exactly our relativistically correct mass relation! 



VI. DIFFERENTIAL CROSS-SECTIONS 



So far, we have considered only total production cross-sections. For calculating differential (instead of total) cross- 
sections in Method I, one must calculate 

der(cc(q)) 
dX 

where A is a kinematic parameter such as the rapidity y, the Mandelstam variable t, or the transverse momentum 
squared "p\. The matching condition is 
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da(cc(q)) _ G n {m c ,-P,x) 



dx m? q) \o) (6.1) 



m c 



and the factorization formula is 

da(H) G n (m c ,P,X) , l/n „. > , . 

Let us contrast this situation to that in Method II. There, the matching condition allows one to calculate the "effective 
Feynman amplitude squared" for charmonium production, and differential cross-sections are obtained in a straight- 
forward manner using that object. There is no need, in Method II, to decide at the level of the matching which sort 
of (differential) cross-section is ultimately desired. 

VII. CONCLUSION 

We have pointed out that there exist two different methods for calculating factorization formulas for charmonium 
production, one suggested by Ref . fl , the other by Ref. || . The results of the two methods can be reconciled using 
the relations 



2M h (O h ) = (O h ) 
2M H (V H ) = (V H ) 



and 



It is interesting to ask: "which method is preferable?" and "Have we truly reconcile d th e two methods?" Indeed, 



one might object to the step in which the equations of motion are applied to obtain Eq. 5.8, since the matrix clement 



of the operator ip' (D 2 %) is ultraviolet divergent and requires a subtraction that is proportional to ip'x 

(V r >) = m c (M r ,-2m c + C) (0") 
where C is a subtraction constant that depends on the renormalization scheme. 
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APPENDIX A: NRQCD CONVENTIONS FOR METHOD I 



In performing the matching procedure, we must evaluate the right-hand side of the matching condition, which is 
written in terms of objects of the NRQCD effective theory. Therefore we require a set of definitions and conventions 
for calculating in this theory. 

Below is presented a practical set of conventions for performing calculations in Method I. These definitions are 
compatible with Ref. which employs the standard nonrelativistic normalization of the quark states. In the heavy 
quark state |c(p, a, i)), a represents spin and i represents color. The states and annihilation-creation operators are 
defined via 



(c(p,a,i)\c(q,T,j)) = (2tt) 3 <5(p - q)5 CTT % 
{a(p,a,i),a t (q,T,j)} = (2tt) 3 <5(p - q)<5 CTT % 

aHp,°,i)\Q) = |c(p,cr, 0) 
a(p,<7,i)|c(q,r, j)) = (2tt) 3 (5(p - q)^ T ^|0) 



The nonrelativistic heavy quark field operators are 







f d 3 k 
(2^ a(k ' ff ' 






=VJ 


j- . Q a' (k, a 
(2?r) 3 v 














"?7 


f (2^' (k '"' 





ik-x 



The free-quark NRQCD operators are 

of q) CS )=xH\cc)(cc\^ X 



3 c c(«)/l 



CS ) 



X t i/'|cc)(cc|V' t (^) X + li.c 



where 



I C (<1> a > *) c (-<ii T > j)) ( c (<ii °"> *)c(-q, t, j)\ 



(Al) 
(A2) 
(A3) 
(A4) 

(A5) 
(A6) 
(A7) 
(A8) 

(A9) 
(A10) 
(All) 

(A12) 



It must be stressed that these operators are intended for use in the calculation of unpolarized production rates, i.e. 
production ra tes in which the a ngular momentum of the final quarkonium state is not specified. 

Using Eqs. Al through A12, one can obtain expressions for the free-quark NRQCD matrix elements, as required 
by the right-hands side of the matching condition. One obtains 

(0|O c "^( 1 5 )|0) = (0\x^Acc)(cc\^x\0) = N^CCtCH* 



xH\cc){cc\ip f (^y) x + h.c. 



\o) = n c J2&Cr£c q 2 



(A13) 
(A14) 



These fr ee-q uark produ ctio n ma trix elements appear in th e rig ht-han d si de of the matching condition given in Eq. 2.1 



and Eq. 2.£. Eqs. A13 and A14 serve in passing from Eq. 2.8 to Eq. 2.9 



The non-perturbative (analytically incalculable) NRQCD matrix elements appearing in the factorization formulas 



<0|0"(^o)|0) = <0| X ty£ \ v + s)( v + s\^ x \o) 



(A15) 



and similarly for (0[P r '( iSb)|0), where S indexes all soft additions to an rj particle, such as one pion, two pions, seven 
pions and a photon, etc., with total energy less than the NRQCD cut-off. 
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APPENDIX B: NRQCD CONVENTIONS FOR METHOD II 



Below is presented a practical set of conventions for performing calculations in Method II. These definitions are 
compatible with the formulations in Ref. |§, which employs relativistic normalization of the quark states. 



<(p, a, i)\c(q, r,i)» = 2E q (2ir) 3 S(p - q) 
{o(p, <J, i),a\bfq, r,j)} = 2E q (2ir) 3 5(p - q) 
a f (p, cr,i) |0» = \c(p,a,i)} 
a(p,a,i)\c(q,T,j)) = (2ir) 3 2E q S(p - q)|0) 



The nonrelativistic heavy quark field operators are 



d 3 k 



(2tt)3 % /2^: 
d 3 k 



6+(k,<7,i)Cae ifc - !B 



(271)3^2^ 
The free-quark NRQCD operators are 



f — ik-a 



xVlc^Xccl^J — ) x + h.c. 



where 



X! I C (<1> °"> ^M -1 !) T ' J'))(c(q, 0", i)c(-q, t, j) | 



(Bl) 
(B2) 
(B3) 
(B4) 



(B5) 
(B6) 
(B7) 
(B8) 

(B9) 
(BIO) 
(Bll) 

(B12) 



Using Eqs. Bl through B12 , one can obtain expressions for the free-quark NRQCD matrix elements, as required by 
the right-hands side of the matching condition. One obtains 



(0|O-M(%)|0» = (0\ X ^\cc}(cc\^x\0} = ^E 2 q N c J2C(rCHa 



(o\r ccM CSo 



xH\cc){cc\ij f x + h.c. 



= AElN c ££CrC& q 2 



(B13) 

(B14) 
(B15) 



These free-quark NRQCD matrix elements appear in the right-hand side of the matching condition given in Eq. B.2. 



The non-pertur 
are 



Eqs. p313| and [B14 serve in passing from Eq. |3.5| to Eq. |3.6[ 

native (analytically incalculable) NRQCD matrix elements appearing in the factorization formulas 



(o\o*>es )\o) = (0|xVE to + S K* + s\^x\o) 

s 

and similarly for (0|7" ? ( 1 5 )|0). 



(B16) 
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